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Based on a study of the analytic running coupling obtained from the standard perturbation theory 
results up to four-loop order, the QCD “synthetic” running coupling a syn is built. In so doing the 
perturbative time-like discontinuity is preserved and nonperturbative contributions not only remove 
the nonphysical singularities of the perturbation theory in the infrared region but also decrease 
rapidly in the ultraviolet region. In the framework of the approach, on the one hand, the running 
coupling is enhanced at zero and, on the other hand, the dynamical gluon mass m g arises. Fixing 
the parameter which characterize the infrared enhancement corresponding to the string tension 
a and normalization, say, at M T completely define the synthetic running coupling. In this case 
the dynamical gluon mass appears to be fixed and the higher loop stabilization property of m g is 
observed. For a = (0.42 GeV) 2 and a syn (M 2 ) = 0.33 ±0.01 it is obtained that m g = 530 ±80 MeV. 

PACS numbers: 12.38.-t, 12.38.Aw. 11.15.Tk 


I. INTRODUCTION 

The paper is devoted to description of the running cou¬ 
pling (the invariant charge) of QCD in which on the basis 
of the perturbative study up to four-loop order the at¬ 
tempt has been undertaken to include the nonperturba¬ 
tive contributions in such a way that nonphysical singu¬ 
larities of the perturbation theory in the infrared region 
be cancelled and the essentials of the QCD dynamics for 
all energy scales accumulated in its framework. 

Doing that it seems notably attractive, on the theoret¬ 
ical grounds of the analytic running coupling of QCD, to 
build the running coupling including additional nonper¬ 
turbative terms. We will consider a variation of such 
couplings, the “synthetic” running coupling a syn (Q 2 ), 
which on the practical grounds of the perturbation the¬ 
ory successfully used for description of the region of large 
Q 2 (short distances physics), contains the nonperturba¬ 
tive terms, determining the main features of QCD in low 
Q 2 region (long distances physics) without the dramatic 
changes of these qualitatively different regimes. 

The analytic approach in QFT was formulated late in 
the 50s in Refs.[l|andil for QED and other theories. For 
QCD the analytic approach was applied in 111 Indi¬ 
cate Refs. ill si ma mm in which the analytic 
approach and its applications to QCD are considering. 
The analyticity requirement, which follows from the gen¬ 
eral principles of QFT, enables one to resolve the dif¬ 
ficulties connected with nonphysical singularities of the 
perturbation theory in the infrared region. In the “an¬ 
alytically improved” running coupling these singularities 
are cancelled out by the nonperturbative contributions. 
In the ultraviolet region the nonperturbative contribu¬ 
tions rapidly decrease and the perturbative contributions 
are decisive. Nonetheless, the behavior of the nonpertur¬ 
bative contributions in the ultraviolet region originating 
from the procedure is of a considerable interest. Their 
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behavior appears to be important in the construction of 
the synthetic running coupling. 

Whereas for the one-loop case a separation of the ana¬ 
lytic running coupling into the perturbative and the non¬ 
perturbative components and behavior of the nonpertur¬ 
bative component are obvious, for the multi-loop analytic 
running coupling this is not the case. For the two-loop 
case such separation in an explicit form and study of 
the nonperturbative component of the analytic running 
coupling was rnadejn Ref. HI for the three-loop case it 
was made in l^dj, 0, and for the four-loop case it was 
made in |16l 1 1 71 . The nonperturbative contributions in 
the analytic running coupling were extracted explicitly 
and their expansion in powers of A 2 /Q 2 was obtained. 
The effective method of the precise calculation of the an¬ 
alytic running coupling was developed on a basis of this 
expansion. 

It is known 3;] that the analytic running coupling is 
finite at zero. In Refs. 0| and Il7l it was shown that the 
finiteness of the analytic coupling at zero was a conse¬ 
quence of the asymptotic freedom property of the initial 
perturbation theory, aan loop (0) = 47 t/ 6 0 ~ 1.396. The 
running coupling is finite at zero also in the case of “freez¬ 
ing” of the interaction ,18|. The loss of Q 2 -dependence 
at low Q 2 takes place for a number of different defini¬ 
tions of the running coupling 0. However, with such 
infrared behavior of the running coupling the descrip¬ 
tion of the confinement and the dynamical breaking of 
the chiral symmetry is not immediate. The behavior of 
the running coupling corresponding to the infrared en¬ 
hancement of the interaction seems to be appropriate for 
description of these most important properties of QCD. 
In general the running coupling constant is ambiguously 
determined. In the perturbation theory it depends on the 
renormalization scheme choice. Besides it can depend on 
the nonperturbative contributions in an observable used 
for definition of the running coupling. In the so-called 
nonperturbative V-scheme for the running coupling the 
singular behavior ay ~ l/Q 2 at Q 2 —* 0 corresponds to 
the linear confining quark-antiquark static potential at 
large distances with the universal string tension param- 
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eter a. Point at Ref. OS where the behavior of the run¬ 
ning coupling of the form ~ 1/Q 2 for all Q 2 takes place 
for a perturbative treatment of theories with permanent 
confinement. In QCD the static potential is defined in 
an explicit gauge invariant form through the vacuum ex¬ 
pectation value of the Wilson loop. Its calculation on 
the lattice in quenched approximation at large distances 
r completely corresponds to the string picture of heavy 
quarks interaction |2li 12211 with confining term err and 
a q/r correction with a coefficient 7 as predicted by the 
bosonic string theory. The synthetic running coupling 
under consideration also belongs to the singular type cou¬ 
plings, a syn ~ l/Q 2 at Q 2 —> 0, and yet it has its own 
motivation related to the study of the nonperturbative 
contributions at Q 2 —> 00 . The synthetic running cou¬ 
pling appears practically to coincide 12 .'ll with the initial 
perturbation theory running coupling (we use a s in the 
MS-scheme) in the region of its applicability but is well 
defined for all Q 2 > 0. 

The main methods of nonperturbative study of the 
Green’s functions in QCD and the running coupling 
which can be built out of these functions are solving the 
Dyson-Schwinger (DS) equations and lattice calculations 
of the functional integrals. In Ref. [3 a summary of the 
results of such recent studies is given, which can be sup¬ 
plemented with papers j25|, [2(| where analytic methods 
were used and |27| with the lattice stimulation results. 
The variety of the results for the behavior of a s (Q 2 ) in 
the infrared region is connected in particular with dif¬ 
ferent truncation methods applied to close the DS equa¬ 
tions. Besides, solving the closed integral equations (or 
systems of equations) requires, as a rule, the simplifying 
assumptions quite often breaking the gauge symmetry 
and the pure technical approximations. It is not sur¬ 
prising, than, that the results of the infrared behavior 
study of a s differ greatly from one another and should 
not be compared literally. Let us note the review pa¬ 
pers [HEED] of the investigations on the IR behavior 
of the Green’s functions, the running coupling in QCD, 
and their applications in the hadron physics. 

The possibility of the singular behavior a s ~ 1/Q 2 at 
Q 2 —> 0 which we consider has been studied in a number 
of papers. In particular, in Ref. |jy] the infrared behav¬ 
ior of the gluon Green’s functions was studied by the 
analytical calculations of the corresponding Feynman in¬ 
tegrals of the DS equation for the gluon propagator in 
the ghost-free axial gauge where the running coupling 
was defined by the full gluon propagator. It was shown 
that the singular behavior of the gluon propagator of the 
form D(Q) ~ 1/(Q 2 ) 2 at Q 2 —■> 0 was possible but it is 
essential to give up the commonly used approximation of 
the three-gluon vertex function by its longitudinal part 
and to take into account the transverse part of the three- 
gluon vertex function of a definite form. 

This paper is organized as follows. In Section QH the 
one-loop synthetic running coupling model of QCD for 
all Q 2 is considered. In Section ITTT1 we study the analytic 
running coupling for the standard perturbation theory 


approximations up to four-loop order and its separation 
into the perturbative and the nonperturbative compo¬ 
nents. In Section ITvI the synthetic running coupling with 
the nonperturbative contributions suppressed at large Q 2 
is build on the basis of the analytic running coupling. 
Setting the parameters of the synthetic running coupling 
is made. In the concluding section the main results are 
summarized and some remarks made. 


II. ONE-LOOP SYNTHETIC RUNNING 
COUPLING OF QCD 


Let us consider the following additive modification of 
the one-loop running coupling of QCD by means of the 
nonperturbative pole type terms 1411 


1 A 2 

ln(Q 2 /A 2 ) + A 2 - Q 2 
cA 2 (1 — c)A 2 

+ ~Q 2 + Q 2 + m 2 J ’ 

where the gluon mass parameter 

A 



( 1 ) 

( 2 ) 


Here Q 2 is the Euclidean momentum squared, constant 
60 = 11 — 2n^/3 (rif is the number of active quark fla¬ 
vors), A is the dimensional parameter of the one-loop 
model 0 , c is the dimensionless parameter of this model 
(it is convenient to introduce the dimensional parame¬ 
ter Ai = t/cA), c £ (l,+oo). The parameter A can 
be fixed, for example, by the normalization condition at 
large Q 2 , whereas the parameter c of the model, as can be 
seen further, describes the relation between the parame¬ 
ter Aqcd and the string tension parameter a of the string 
models [50|. It stands to reason that for the realistic def¬ 
inition of the parameters Aqcd, c and their connection 
it is necessary to go out of the one-loop approximation. 

The first term of Eq. m is the solution of the renor¬ 
malization group equation for the QCD running coupling 

a s (Q 2 ) 


Q 


,da s (Q 2 ) 
dQ 2 


= 0{a s ) 


( 3 ) 


in the one-loop approximation, j3(a a ) — —fi^a 2 , /3 q = 
bo/4n. Introducing the renormalization invariant param¬ 
eter A (the integration constant of the differential equa¬ 
tion) we obtain 


a 




47T 1 

Hq ln(Q 2 /A 2 ) ’ 


( 4 ) 


where a nonphysical singularity (the Landau pole) at 
Q 2 = A 2 is present. The vanishing of expression © 
for Q 2 —> 00 corresponds to the remarkable property 
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of asymptotic freedom 39] of non-Abelian gauge theo¬ 
ries, while the growth of a s (to some critical value or to 
infinity) with decreasing Q 2 can be connected with the 
confinement problem. 

The pole terms in Eq. 0 are nonperturbative, A 2 ~ 
p 2 exp {— 4-rr/boCXs(p 2 )}- The sum of the first two terms 
in Eq. © is an analytic function in the complex Q 2 -plane 
with a cut from 0 to —oo, 


v(l) 


m = 


47T 


ln(Q 2 /A 2 ) 


A 2 ' 
A 2 -Q 2 


(5) 


This function can be presented by the dispersion relation 


ain(< 3 2 ) = J 

0 


da p^\a) 
a + Q 2 


( 6 ) 


where the function p W (cr) is called the one-loop spectral 
density 


p (1) (a) 


47T 1 

bo ln 2 (cr/A 2 ) + 7 r 2 


(7) 


The last equation can be obtained by the analytic contin¬ 
uation of a s into the Minkowski space Q 2 —> — cr — *0 and 
calculation of the imaginary part p^\a) = (—a — 

*0). For real Q 2 > 0 function © is positive monotone de¬ 
creasing function with maximum at zero ain ( 0 ) = 47t/&o- 
The second nonperturbative term in Eq. © does not con¬ 
tribute to the imaginary part of ain (Q 2 ) in going to the 
Minkowski space, so that p^\a) = ^3ain (—a — iO). 
The synthetic running coupling © can also be presented 
in the form of the dispersion relation 



da )0syn(<7) 
a + Q 2 


( 8 ) 


The function piyl(a ) will be called the one-loop spectral 
density for the synthetic running coupling. It contains 
additional terms in the form of the delta functions, 

piyn(o-) = /5 (1 V) + ^ [cA 2 6(a) + (1 - c)A 2 <5(ct - m 2 )] . 

(9) 

Introducing two pole terms at Q 2 = 0 and Q 2 = —m 2 < 0 
does not change the analyticity domain of the analytic 
running coupling ©. In Eq. © for the spectral density 
the additional terms in the form of two ^-functions local¬ 
ized at a = 0 and a = m 2 > 0 emerged (in comparison 
with expression © of the perturbative origin). 

Let us bring equation © for aiyn (Q 2 ) to the explicit 
renormalization invariant form. It can be done without 
solving the differential renormalization group equations. 
Writing the normalization condition for c4yn(Q 2 ) we ob¬ 
tain an equation for the required dependence of the pa¬ 
rameter A 2 on the values asyli(/i 2 ) and p 2 of the form 

m , 2 \ 47 T f 1 A 2 cA 2 

synlM ; bo [ln(/x 2 /A 2 ) A 2 - p 2 p 2 


, (1 ~ c)A 2 ~ 

p 2 + TO 2 

From dimensional considerations 


( 10 ) 


A 2 = p 2 exp{-</? (a(/i 2 ))}, (11) 

where a(p 2 ) = ( 6 0 / 47 r)aiyn(/x 2 ). Then for ip (a) we have 
a transcendental equation 

a=-^ + --—-pr + ce _v(a) — - ( . C ~^ 2 ( 12 ) 

ip (a) 1 — 1 + (c — l)e^( a l 

The function ip(a) has the behavior ip (a) ~ 1/a —■> +00 as 
a —> +0 for all values of c. This behavior corresponds to 
the perturbative region. The behavior of this solution 
at a —> +00 is ip (a) ~ — ln(a/c) —> —00. The beta 
function /3 S y n (<a S yn) for the synthetic running coupling 
can be found by the equation which is analogous to the 
equation 

Q 2 ~ |4g) =AyB(a , y „). (13) 

Differentiating the running coupling © with the use of 
equations o and dm we obtain 


4:7 r 

/^syn(ctsyn) = 


—a + 


1 

<p(a) 


1 1 

p 2 (a) (1 — e v f a l ) 2 


(c- 1) 2 1 

(l + (c- l)e ^(“)) 2 J 


a=bo«ayn/4ir 


(14) 


Therefore, using the behavior of function p(a) as a —> 0 
and 00 , we find the asymptotic behavior 


/^syn(cfcsyn) 


bo 2 

4^ syn 


+ °( a syn)i 


*syn 


0 , 


(15) 


/^syn(o^syn) — C^syn c(c 2)-t-(?(l), Ctsyn » OO. (16) 
DO 

Doing the corresponding expansions we make sure that 
the singularity of the /J-function m at ip —> 0 is seeming. 
We also make sure that for all a syn > 0 the function 
/3syn(a S yn) is negatively defined. 

Let us write the last three terms of the synthetic run¬ 
ning coupling ©, taking into account ©, in the form 


<y P ^ ( 1 ) (Q 2 ) 


47t A 2 

h [A 2 ~Q 2 


cA 2 

~Q 2 


(1 — c)A 2 ’ 
Q 2 + m 2 _ 


4 - 7 T cA 6 

5o" Q 2 (A 2 — Q 2 )(A 2 + (c — 1)Q 2 ) ’ 


(17) 


For a —> +0, according to equation m the function 
(j>{a) ~ 1/a, thus for A 2 from equation 1111 ) we obtain 
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A 2 ~ /x 2 exp (—1/a) at a syn (/r 2 ) = (47r/5 0 )a —*• +0 so 
expression ra must be considered as a nonperturbative 
component of the synthetic running coupling. The be¬ 
havior of the nonperturbative “tail” at large Q 2 is the 
following: 


a 


npt(l) 

syn 


(Q 2 ) 



c A 6 ' 

7^1 Q*_ 


+ o 



(18) 


As seen from Eq. ca, the nonperturbative contributions 
of the synthetic running coupling decreases at large Q 2 
substantively faster than that of the analytic running 
coupling ©• 

Hence the one-loop synthetic running coupling of QCD 
has the following interesting properties: 

(i) By the construction, as a function of Q 2 , it has an 
analytic structure corresponding to the causality; 
that is, it is a holomorphic function in the com¬ 
plex Q 2 -plane with a cut along the negative real 
semiaxis. 


(ii) As a function of its value a syn (/i 2 ) at the normaliza¬ 
tion point /i 2 , it has an essential singularity at the 
origin; the asymptotic expansion of its nonpertur¬ 
bative part in a syn (fi 2 ) for a syn (/.i 2 ) —> +0 is equal 
to zero, which ensures conformity to the initial per¬ 
turbation theory. 

(iii) In the ultraviolet region, it coincides with usual re¬ 
sult of the perturbation theory (with renormaliza¬ 
tion invariance taken into account) apart from fast 
decreasing power terms. The nonperturbative com¬ 
ponent behaves as ~ 1/(Q 2 ) 3 when Q 2 —» oo. 

(iv) In the infrared region the synthetic running cou¬ 
pling does not have nonphysical singularities of the 
perturbation theory. There is the mass term and 
the singular at the origin term which can be re¬ 
sponsible for the confinement of quarks. 


For small Q 2 ( Q 2 < A 2 ) we have the expansion 


4-7T A 2 


s'W'^EMi-nis 


n —0 




A 2 


( 21 ) 

The case c = 2 is the particular one from the symmetry 
considerations, 


a 


npt (1) 
syn 


(Q 2 ) 


47T 

b 0 


‘ A 2 
A 2 -Q 2 


2A 2 

~Q 2 


A 2 ‘ 
A 2 + Q 2 


_ Hr 2A 2 _ A* _ 

bo Q 2 A 4 - Q 4 ’ 1 j 

for which the nonperturbative component is the odd func¬ 
tion of Q 2 , 


<yn (1) (-Q 2 ) = -«syn (1) (Q 2 )- (23) 

Here, in the ultraviolet expansion as well as in the in¬ 
frared expansion there are no terms of the even powers 
of Q 2 . In particular, the infrared expansion does not have 
the Coulomb’s mode. 

Indicate for completeness two boundary cases for the 
values c G (l,+oo) considered. The first case is c = 1, 
for which 


a 


npt(l) 

syn 


(Q 2 ) 


47T A 4 
fro" Q 2 (A 2 — Q 2 ) ’ 


and the nonperturbative component decreases at infinity 
not so fast as in expression ®. The second one is c = oo 
for which 


«sy P ^ (1 ) (Q 2 ) 


4tt A 6 
fro" Q 4 (A 2 — Q 2 ) ’ 


and the singularity in the infrared region is stronger, 
Q!syn(Q 2 ) ~ 1/(Q 2 ) 2 , Q 2 - 0. 


As it will be evident from the subsequent considera¬ 
tions, all these properties are valid for the two-, three- 
and four-loop synthetic running coupling. It is signifi¬ 
cant that for the one-loop case we have not only rep¬ 
resentation ©, but for one thing, the nonperturbative 
contributions are extracted from the synthetic running 
coupling in the explicit form 

4n(Q 2 ) = a pt(1) (Q 2 ) + < n t(1) (Q 2 ) (19) 

and for another, for the nonperturbative contributions, a 
simple formula is on hand. For large Q 2 (Q 2 > A 2 ) this 
contributions can be represented as a series 


III. MULTI-LOOP ANALYTIC RUNNING 
COUPLING OF QCD 


For the multi-loop case the renormalization group 
equation © for the QCD running coupling a s (Q 2 ) is 
of the form 


Q 2 


da s (Q 2 ) 

dQ 2 


P(<*s) 


= -Poa 2 s ~ Pia 3 s - fca* - fi 3 a 5 s + 0(a 6 s ). (24) 


a 


npt (1) 
syn 


(Q 2 ) 


47T 

'h 


A 2 

Q 2 


3 oo / 


n—0 


A 2 V 

vq 2 ; 


( 20 ) 


The coefficients /3o, Pi do not depend on the renormal¬ 
ization scheme choice, whereas the next coefficients do 
depend on this choice. For the numerical calculations we 
use its values within the MS- scheme. 
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Let us write down the solution of equation (12411 for 
a s {Q 2 ) at L = ln(Q 2 /A 2 ) —> oo in the form of the stan¬ 
dard expansion in inverse powers of the logarithms 




PoL 


(3 2 L 


of the logarithmic function). In Refs. IBE3 up to four- 
loop order the separation of the analytic running coupling 
into the perturbative and nonperturbative components 
was obtained, 

a an (Q 2 )=a pt (Q 2 )+aZ\Q 2 ). (30) 



Pi 

In 2 T 

In L - 

1 1 

PoPi 



' PtL 2 



ft J 



Pi 

In 3 L 

-®In= 

L — 

^2 — 

SPoPl 

j In L 


PpL 3 


2 



PI 



1 

+ 2 _ 

PM 

2p 3 \ 

+ o 

(i 

)}• 


(25) 

The sum of the terms of 

Eq. ED 

up to 

1/L n 

order 


(n=l,2,3,4) will be referred to further on by the n-loop 
perturbative component of the running coupling and de¬ 
noted as a pt (Q 2 ). It can be written in the form 

A.'TT 

a pt (Q 2 ) = —a(x), (26) 

bo 


In Eq. m we take a pt (Q 2 ) as the initial (in our case the 
standard) solution of the renormalization group equation 
up to four-loop order and aO pt (Q 2 ) appears to arise ad¬ 
ditionally as a result of the procedure. The following 
expansion was obtained in the power series 



where the coefficients c n were defined by the beta func¬ 
tion coefficients. Note the important properties of ex¬ 
pansion GD such as the higher loop stability of the co¬ 
efficient of the leading term and the slow increase of the 
coefficients c n with number n. 


IV. MULTI-LOOP SYNTHETIC RUNNING 
COUPLING OF QCD 


a(x) 


1 ^ln(lna;) In 2 (In a;) 
lux In 2 a: [ In 3 a; 


ln(lna;) 
In 3 x 


K 

In 3 x 


-b 3 


In 3 (In a:) 
In 4 x 


5 In 2 (In a;) Nn(lnad 

2-nyr + (3 ' t+1, i!y J 


X 


Here x = Q 2 /A 2 and the coefficients 


In 4 x 
(27) 


The one-loop synthetic running coupling can be nat¬ 
urally extended to the multi-loop cases. Thus mod¬ 
ify the analytic running coupling introducing two addi¬ 
tional nonperturbative terms, the singular at zero term 
of the form ~ 1/Q 2 and the mass term of the form 
~ 1/(Q 2 + TO 2 ). As a result, we come to the expression 


i{Q 2 ) — «an(Q 2 ) + 


47T 

bo 


IQ 2 


dA 2 ' 
Q 2 + m 2 \ ’ 


(32) 


b 


Pi 


K = — 1 + 


ftp Pi 

PI 


- _ 1 PlPz 

K 2 2/3? ' 


(28) 


The values of b 0 , b , k, k depend on rif. Within the 
standard picture of matching the solutions at heavy 
quark thresholds the parameter A becomes dependent 
on n, f [rT| . 

The analytic running coupling a an (Q 2 ) = 
(4n/bo)a a . n (x) is defined through the dispersion re¬ 
lation 


OO 

a an (x) = - [ p(cr), (29) 

7T J X + <7 
0 

where the spectral density p(a) = 3a an (—cr — *0). The 
analytic approach suggests[52j that 3a an (— o — *0) = 
3a(—cr — *0). As a result from function a(x) of the 
form (ED with nonphysical singularities on the positive 
real semiaxis of the complex plane x = Q 2 /A 2 we come 
to the function a an (x) of the form (12911 , which is a single¬ 
valued analytic function in the complex plane x with a 
cut from 0 to — oo (with a standard definition of the cuts 


containing, besides A, three subsidiary parameters: c, d 
and m g ( m g = toaA), which are assumed to be nonzero. 
We will define this parameters in the following way. At 
large Q 2 ( Q 2 > A 2 ), using expansion ED , we can find 


a S yn(0 2 ) = aP t (0 2 ) + 


4-7T 


A 2 

(ci + c-d)^ 


+(c 2 + dm\) 



(c 3 - dm\) 



3" 



Demand the nonperturbative contributions to be mini¬ 
mal at large Q 2 , i.e. the terms of the form ~ 1/Q 2 , 
~ l/(<5 2 ) 2 to be absent in Eq. (13311 [ S3 | . Then two of 
three parameters are fixed by the following equations 

d=c + ci, m\ = -c 2 /(c + ci). (34) 
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The parameter Ai = y/cA will be considered as fixed. 
The coefficients c n < 0, therefore the absence of tacliion 
condition is A < A\/y/—c\. With a given number of loops 
the free parameter of the synthetic running coupling is 
only one parameter A. Then taking into account Eq. m 
we have 


a syn from the correspondence of the singular at zero term 
in Eq. © to the infrared behavior itHHll of the running 
coupling ay Then 

b = tlAf, A, = cA 2 . (40) 

Z Do 


a syn (Q 2 ) = a pt (Q 2 ) + — 
o o 


c 2 A 2 1 
A 2 + Cl A 2 _ 



+0((Q 2 ) 4 ). (35) 

As seen from Eq. m the leading power nonperturbative 
term decreases rapidly at Q 2 —> oo and in the absence of 
the tachion is negative l54j . For the gluon mass param¬ 
eter m g which we call the dynamical gluon mass from 
Eqs. (13411 it follows that 


Therefore, if the string tension parameter is given, the 
parameter Ai can be specified by Eq. m . Then with 
a ~ 0.42 GeV, b 0 = 9 we obtain^ Ai ~ 435 MeV. The 
parameter A (as well as the parameter c) can be fixed by 
the normalization condition, and the synthetic running 
coupling will be fixed completely. 

Consider the dependence of the dynamical gluon mass 
on A for different number of loops of the initial perturba¬ 
tion theory approximation. In Fig.^the dynamical gluon 
mass m g { A) is shown for one — four-loop cases. Up to 
400 MeV the curves do not diverge too much and at A = 
375 MeV m g ~ 0.6 GeV. We normalize the running cou¬ 


m.g = A A 


I -c 2 A 2 
A 2 + dA 2 ' 


(36) 


Let us turn to the interpretation of the parameter Ai 
which describes the value of the singular term. As it has 
already been noted, the behavior of the running coupling 
ay ~ 1/Q 2 at Q 2 —> 0 corresponds to the linear confining 
quark-antiquark static potential in quenched QCD. We 
set up a correspondence of the potential and the running 
coup ling using the equation (to compare with Refs. |4l| 
and 13611 




4nay(q 2 ) 


(37) 


77.-3 


where ay{q 2 ) is defined as an effective charge which 
is the renormalization scheme independent and gauge in¬ 
variant quantity. The color factor corresponds to the 
SU(N C ) group, N c = 3. Let us assume that in the in¬ 
frared region 


a v {q 2 ) 


3 a 

2 ?’ 


0 . 


(38) 


Then the integral over three dimensional momentum 
space in Eq. m formally diverges at the origin. We de¬ 
fine this integral introducing the dimensional regulariza¬ 
tion. After integration over the n-dimensional Euclidean 
momentum space we put n = 3. Then as the divergences 
do not occur, the transition to the n-dimensional integra¬ 
tion provides not only the regularization but the defini¬ 
tion of the divergent integral for the case n = 3. As a re¬ 
sult, for the infrared behavior of the effective charge 01 
the behavior of the potential at large distances is as fol¬ 
lows: 


V(r) ~ or, r —> oo, (39) 

where a = a 2 is the string tension parameter. Let us 
define the parameter Ai of the synthetic running coupling 


TABLE I: The parameters A pt (MeV), A an (MeV), A syn 
(MeV), the dynamical gluon mass m g (MeV) and the param¬ 
eters c, d on the number of loops. The number of active quark 
flavors rif = 3, Ai = 435 MeV. The normalization condition 
is a(M 2 ) = 0.33, M T = 1.777 GeV. 



l-loop 

2-loop 

3-loop 

4-loop 

Apt 

214.25 

395.10 

364.19 

357.32 

Aan 

254.51 

636.02 

523.86 

535.54 

Asyn 

214.26 

397.10 

365.26 

358.70 

m g 

121.14 

648.09 

461.66 

526.57 

c 

4.1282 

1.2018 

1.4204 

1.4728 

d 

3.1282 

0.5358 

0.7706 

0.7447 


plings Qs yn (O ' 2 ), a an (<3 2 ) and a pt (Q 2 ) at the r-lepton 
mass by QEB| a{M 2 ) = 0.33, M T = 1.777 GeV. For 
this normalization condition the values of the parameters 
A, the dynamical gluon mass m g and the parameters c, 
d, are given in Table 1 Point to two things. The pa¬ 
rameters A syn and A pt are close in value whereas the 
parameters A an are considerably larger. This is a conse¬ 
quence of conditions (El which give the fast decrease of 
the nonperturbative terms of a syn (Q 2 ) at large Q 2 . For 
all quantities considered the one-loop case turns out to be 
exceptional, and then the stabilization is observed with 
the number of loops of the initial perturbation theory 
approximation. 


V. CONCLUSIONS 

In the construction of the QCD analytic running cou¬ 
pling (1291) the nonphysical singularities of the perturba¬ 
tion theory in the infrared region disappear and in the 
ultraviolet region the nonperturbative power corrections 
arise decreasing rapidly at large Q 2 , in comparison with 
the main perturbative component. However, when con¬ 
sidering the nonperturbative quantities it may happens 
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FIG. 1: The gluon mass parameter m g as a function of A 
approximation. 

that the decrease of the nonperturbative contributions 
is not fast enough for the consistent definition of these 
quantities. In the synthetic running coupling it is pro¬ 
posed to provide the highest possible suppression of the 
nonperturbative contributions at large Q 2 by means of 
a minimal number of the pole type terms. The parame¬ 
ters characterizing the additional nonperturbative terms 
have a clear physical meaning and take the reasonable 
values. The running coupling 0 built on the basis of 
the analytic running coupling is called the one-loop 
synthetic running coupling because it contain the param¬ 
eters related to the ultraviolet region as well as to the 
infrared region. We introduce the singular at zero term 
of the form ~ 1 /Q 2 which can correspond to the lin¬ 
ear quark confinement and the mass term of the form 
~ 1 /{Q 2 + to 2 ) with the parameter m g corresponding to 
the non-vanishing dynamical gluon mass. We impose the 
condition of the fastest decrease of the nonperturbative 
component at large Q 2 and receive the running coupling 
model of form ®. The model asyn(Q 2 ) has two indepen¬ 
dent parameters, the dimensional parameter A and the 
dimensionless parameter c which defines the value of the 
singular term. In Section 2 the one-loop model of the syn¬ 
thetic running coupling and its nonperturbative compo¬ 
nent properties for c € [l,+oo) are considered. A study 
of the multi-loop analytic running coupling and its non¬ 
perturbative component provide a possibility of natural 
generalization of the synthetic running coupling model 


different number of loops of the initial perturbation theory 


to the multi-loop cases. The multi-loop synthetic run¬ 
ning coupling lTT21l as the one-loop model is constructed 
by introducing two additional nonperturbative terms of 
the form ~ 1/Q 2 and ~ 1 /(Q 2 + to 2 ). The minimal¬ 
ity principle of the nonperturbative contributions in the 
perturbative region leads to two equations m for the 
introduced nonperturbative parameters. As a result, the 
synthetic running coupling has two independent param¬ 
eters. First, the parameter A which owing to highly fast 
decrease of the nonperturbative contributions at large 
Q 2 practically coincides with the parameter Aqcd m the 
region of application of the perturbative solutions. Sec¬ 
ond, the dimensionless parameter c (or the dimensional 
parameter Ai = y/cA) determining the value of the sin¬ 
gular term. Correlating this parameter responsible for 
the infrared enhancement by Eqs. m with the string 
tension parameter a of the string models, we arrive at 
the dynamical gluon mass as a function of the parameter 
A defined by Eq. (I.'Fill . The corresponding dependencies 
for 1-4-loop cases are given in Fig. Q] with Ai = 435 MeV 
(that is ct 1 / 2 = 0.42 GeV). The normalization completely 
defines the synthetic running coupling and for a(A'I 2 ) = 
0.33 the values of the parameters are shown in Table HI 
The parameter m g in Eq. m for A < A\/y /—ci is real, 
and the synthetic running coupling a syn (Q 2 ) in Eq. 1321) 
is a holomorphic function in the complex plane Q 2 with a 
cut along the real negative semiaxis. Thus the synthetic 
running coupling of QCD a S yn(Q 2 ) has the properties 
outlined in Section 2 for the one-loop synthetic running 


for 









coupling model. 

As seen in Table|I]the parameters of the synthetic run¬ 
ning coupling fixed in such a way show the higher loop 
stabilization. In particular, the dynamical gluon mass 
Trig can be estimated as 400-600 MeV. For a = (0.42 
GeV) 2 and a syn (M 2 ) = 0.32, 0.33, 0.34 it is obtained 
that m g = 453, 527, 613 MeV (for the 4-loop case) [56|. 
Hence the string tension identification of the parameter 
of the synthetic running coupling defining the value of 
the singular term results in the consistent values for the 
other parameters considered. 

According to Eq. (EI51) the nonperturbative contribu¬ 
tions decrease at large Q 2 as ~ 1/(Q 2 ) 3 , which is suffi¬ 
cient for the convergence of the gluon condensate in the 


ultraviolet region |35l |. Thereupon the generalization of 
the gluon condensate studies |H, to the multi-loop 
synthetic running coupling is of much interest. 
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